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PREFACE 


la  June  1991.  Det  1.  4  WW  (aow  50  OSS/WE),  »ked  the  United  Sii#e«  Air  Force  EavironmenuJ 
Technical  Applictaioai  CeoU»'  (USAFETAC).  to  ttady  past  Global  Pofiboatai  Systena  (GPS)  satellite 
anomalies  and  determiae  if  a  statical  correlsatoo  coold  be  fouad  between  commonly  observed  space 
miviroomental  data  and  the  GPS  tatollim  •nomalies.  The  n^nest  (performed  tiader  USAFETAC  project 
number  910824)  was  assigned  to  USAFETACS  Simulation  mid  Techniquet  Bimch  (SYT)  This  i^Kjrt 
documents  the  work  done  on  this  project  and  the  results,  which  are  of  interest  to  the  space  envtronisimiai 
communitv  in  general,  md  in  particulm’  to  those  who  support  GPS 

USAFETAC.'SYT  analysts  compmed  distributions  of  space  environmental  data  to  d«enaine  if  there  were 
differences  between  the  environment  during  nomalies  and  the  general  space  environment  and  calculated 
correlation  coelTicients  for  space  environmentai  variables  We  used  stqiwtte  lincaur  regretsioa  and 
discriminant  analysis  to  determine  which  environmental  elements  had  the  greatest  innuenot  on  the 
occurrence  of  GPS  anomalies,  nd  how  successfully  these  elements  could  disungutsb  cases  in  the  past 
when  GPS  anomalies  occurred 

The  author  wishes  to  thank  Mr  J.  Berg  and  the  other  employees  of  Aerospace  Corporation,  and  Capt 
C.  Larcomb.  SSDA^'E,  for  their  assistance  in  telephone  coosultatioBS  on  GPS  and  in  acquiring  GPS  orbital 
data  as  well  as  the  tracking  models  DGEN  and  PCSOAP.  The  author  also  wishes  to  thank  Capt  G  Deuel, 
National  Geophysical  Data  Center,  for  his  help  in  acqumng  Geostationary  Operational  Environmental 
Satellite  (GOES)  x*ray  and  energetic  proton  data  for  this  study. 
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1.  INTRODUCTION 


1.1  Purf>os«  of  Study.  This  study  undi^ook  to 
fiod  correlations  (if  any)  between 
soiar/geom  luetic  activity  and  the  occurrence  of 
Global  Positioning  System  (GPS)  satellite 
anomatira  (called  "bit  hits"  because  of  their  effects 
of  satellite  software).  When  GPS  satellite 
miomalies  occur,  space  environmental  forecasten 
at  the  AF  Space  Forecast  Center  (AFSFC) 
examine  solm,  geom^netic,  and  siaellite  particle 
data  to  determine  if  the  space  environment  caused 
the  anomaly.  Thoe  has  bemi  no  comprehensive 
study,  however,  to  determine  if  there  is  a  pattern 
between  certain  types  of  activity  in  the  space 
environment  and  the  occurrence  of  GPS  satellite 
anomalies  over  time.  This  study  examines 
commonly  observed  space  environmmital  facton 
and  provides  guidance  on  how  influential  they  are 
in  the  occurrence  of  GPS  bit  hit  anomalies. 

1.2  Limitations.  The  number  of  GPS  anomalies 
reported  during  the  7  years  covered  by  this  shidy 
was  relatively  small,  making  the  occurrence  of  a 
satellite  anomaly  at  any  one  time  a  relatively  rare 
event.  Such  a  small  sample  size  increases  the  risk 
that  a  few  bad  data  points  or  extreme  cases  in  the 
anomaly  data  could  have  a  disproportionate 
influence  on  the  results. 

The  regression  equations  used  to  predict  the 
probability  of  anomaly  events  are  closely  tied  to 
the  database  used  to  derive  than.  When  different 
time  periods  within  the  7  years  of  the  study  were 
used  to  derive  the  equations,  the  results  sometimes 
varied  a  great  deal.  No  independent  dataset  was 
available  to  test  how  well  the  equations  could 
predict  die  probability  that  an  anomaly  would 
occur. 

To  determine  which  environmental  variables  could 
distinguish  between  anomaly  cases  and  those 
hours  without  anomaly  reports,  we  had  to  assume 
an  anomaly  did  not  occur  when  no  report  was  in 
the  database.  Since  indicMions  imply  the  anomaly 
dtdabase  is  incomplete,  our  technique  may  perform 
differently  than  the  skill  indicators  provided. 


1.3  Methodology.  GPS  anomaly  reports  were 
combined  with  ail  3<>iuls^ie  space  envu-onmeaial 
data  to  creMe  databases  These  databases  were 
sorted  into  a  frequency  distribution  for  each 
environmental  vanable  so  that  the  type  of  space 
environment  occurring  dunng  anomalies  could  be 
compared  to  the  general  space  oivironmeni 
Statistical  analysis  was  also  dose  on  the  dai^ases 
to  determine  which  enviroiunenLai  variables  might 
significantly  influence  the  occurrence  of  saieilite 
anomalies 

1.4  RotultS.  Our  results  are  described  in  tbrt% 
stations . 

•  Results  from  comparisons  of  the 
enviroiunental  variable  distributions 

•  Correlation  coefficients  between  the 
environmental  variables  and  the  occurrence  of 
satellite  anomalies, 

•  Results  ftoro  discriminant  analysis  using 
linear  regression  equations  to  classify'  the  data  into 
yes/no  cases,  based  on  the  probability  of  an 
anomaly  occurring. 

Figure  1  contains  regression  equations  for  seven 
groups  of  environmental  variables  discussed  in 
Section  5  of  this  study.  To  classify  die  data  by 
cases  in  which  anomalies  were  likely  or  unlikely 
to  occur,  the  probability  of  anomaly  occurrence 
was  calculated  using  the  results  of  die  r^ression 
equations  in  this  formula; 


Prob  {Yes) 


\ 

exp  {No-Yes)  *  1 


111 


The  probability  of  .5  or  greater  was  used  to  infer 
that  an  anomaly  occurred.  The  equations  for  the 
second  group  of  variables  seem  to  perform  the 
best  of  the  seven  attempts  Further  details  are 
provided  in  Section  5. 
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Qmm  1:  Eriyiteiitiiiiiii  it  of  anomaly  (I98€.9it 
Y««  «-3.24  ■K3.0243*MINSO  +0.0345*F10  -4524*XRAY  +0  OOOOrMfN.2 
•►6.03E-7*MIN95  -0.0049*MIN10  -0.0008*MAX50  +0.000rMAX10  +0.0015*AVG5 
No  s-4.17  +0.0017*MIN50  •►0.0438*F10  -8711*XRAY  +9.16E-6‘MtN.2 
+8,69E-7*MIN95  -0.0017*MIN10  ••■0.0007-MAX50  -O.OOOS'MAXIO  +0.0007*AVG5 

Group  2:  Envtronmottt  it  timt  of  «oofinlv  t19»6.9H 
I  Yo*  =-1.65  +0.0214*MIN50  ■K).0274*F10  -4530*XRAY  -0,0022*MIN10 
•*-0.0007*MAX50-0.0010*MAX10  '►0.0012*AVG5 
No  =-2.71  +0.0004*MIN50  ■f0.0372*F10  -7962*XRAY  -0.0007*M1N10 
+0.0002-MAX50  +6.96E-6*MAX10  -t-0.0014*AVG5 

Group  3:  Envtronmpnt  24-hour»  boforo  onoiwplv  (1986-911 

Y«8  =-3.10  +0.035*F10  -O.OOOrMINI  +0.0052*AP24  +0.000012*M1N.2  -5515‘XRAY 

+4.281 9E-9*AVG30 

No  =-4.11  +0.045*F10 -0.0001*MIN1  +0.0014*AP24  +0.00001  rMrN.2 -8838*XRAY 
+2,981 2E-9*AVG30 

Group  4;  Environmont  48-hour»  boforu  opomotv  (1986-81) 

Yos  =-3.26  +0.034*F10  +0.0006*MIN1  +0.000012*MfN.2  -0.0004*MiN5  -.0004*MIN10 
-0.0002*AVG1  +0.0376*AP 

No  =-4.22  +0.043*F10  -0.0006*M1N1  +0.00001 1*MIN.2  +0.0026*MIN5  -.0034*M1N10 
+0.000rAVG1  +0.0347*AP 

Group  S:  Environiwnt  72.Hour»  bofort  onomHv  (1986-911 

Yos  =-3.24  +0.034*F10  +0.00001 1*MIN.2  +0.0020*MIN50  +0.034*AP  -0.00007*MAX1 

+4.242E-7*MIN95 

No  =-4.24  +0.042*F10  +9.90E-6*MIN.2  -0.0024*MIN50  +0.032*AP  -0.00004*MAX1 
+6.009E-7*MIN95 

Group  6:  Envirowmont  during  soior  min  (1986-881 

Yos  =-3.53  +0.038*F10  +0.0106*MIN50  +7.12E-6*MIN.2  +7.60E-7*MIN95 

-3.2E-9*MAX30  +3.04E-8*AVG30  -2.32E-8*MIN30 

No  =-5.12  +0.082*F10  +0.0007*MIN50  +6  19E-6*MIN.2  +1.07E-7*MIN95 

+5.0E-11*MAX30  +5.29E-9*AVG30  +4.17E-9*MIN30 

Group  7:  Environmont  during  solar  max  (1989-911 

Yos  =-1.30  +0.003*MAX5  -0.0005*MAX50  +0.051 S^AVGSO  +8590*XRAY 

-0.00014*AVG1  +0.0005*AVG10  -0.007*MAX10  -0.0004*MIN10 

No  *-0.02  .0.0005*MAX5  -0.0004-MAX30  +0.0013* AVG50  +4644*XRAY 

+0.00012*AVG1  -0.0006*  AVGIO  +0.0009*MAX10  +O.0002*MINi0 


Ftfare  1.  Regreision  eqaatioss  to  predict  the  probability  of  anomaly  occurrence. 
Environmental  variables  are  defined  in  Figure  2. 


Efivlronniafital  Variabla 

Oita  Type 

Units 

24-Hottr  A, 

AP24 

None 

3-Hour  a. 

AP 

None 

Hourly  30-300  KeV  Electros 

Flux 

MIN30.  AVG30, 
MAX30 

Particles/Steradian/Sec/McV 

Hourly  93-300  KeV  Electron 

Flux 

MIN95 

Paiticies/Steradian/Scc/MeV 

Houriy  .2-2  MeV  Electron  Flux 

MIN.2 

Particlcs/Stcradian/Sec/MeV 

10.7  Cm  Radio  Flux 

FIO 

10  “  Watts/M^/Hcrtz 

X-it^  Flare  Rqiort 

XRAY 

Watts/M’ 

Hourly  >i  MeV  Proton  Flux 

MINI,  AVGl,  MAXI 

Particies/Steradian/Sec/MeV 

Hourly  >5  MeV  Proton  Flux 

MINS,  AVG5,  MAX5 

Particies/Steradian/Sec/MeV 

Houriy  >10  MeV  Proton  Flux 

MINIO,  AVGIO, 
MAXIO 

Paiticles/Steradian/Sec/MeV 

Hourly  >50  MeV  Proton  Flux 

MINSO,  AVG50, 
MAX50 

Particles/Steradian/Sec/MeV 

Figure  2.  Euvironmoital  Variables  Used  in  the  Study.  Data  types  are  defined  below; 


AP24  24-hour  global  geomagnetic  index  (unitless) 

AP  3-hour  global  geomagnetic  index  (unless) 

AVG1  hourly  average  >1  MeV  protons  In  particto^ter/sec/MeV 

AV65  hourly  average  >5  MeV  protons  in  particles/ster/sec/MeV 

AVGIO  hourly  average  >10  MeV  protons  in  partides/ster/sec/MeV 

AVG30  hourly  average  30-300  KeV  eiecbons  in  particles/ster/sec/MeV 

AVG50  hourly  average  >50  MeV  protons  in  particles/ster/sec/MeV 

F10  10.7  cm  solar  radio  flux  in  10  “  watts/m^eite 

MAXI  hourly  max  >1  MeV  protons  in  particles/ster/sec/MeV 

MAX5  hourly  max  >5  MeV  protons  in  partictes/ster/sec/MeV 

MAX10  hourly  max  >10  MeV  protons  in  particies/ster/sec/MeV 

MAX30  hourly  max  30-300  KeV  electrons  in  paiticles/ster/sec/MeV 

MAXSO  hourly  max  >50  MeV  protons  in  particles/ster/sec/MeV 

MIN.2  hourly  min  .2  -  2  MeV  electrons  in  particles/ster/sec/MeV 

MINI  hourly  min  >1  MeV  protons  in  particles/ster/sec/MeV 

MtN5  hourly  min  >5  MeV  protons  in  partides/ster/sec/MeV 

MIN10  hourly  min  >10  MeV  protons  in  particies/ster/sec/MeV 

MIN30  hourly  min  30-3CK}  KeV  electrons  in  particies/ster/sec/MeV 

MIN50  hourly  rrvn  >50  MeV  protons  in  particles/ster/sec/MeV 

Mlh^  hourly  rrnn  95-300  KeV  electrons  in  particles/ster/sec/MeV 
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2.  DATA 


2.1  Sateilit*  Anomaly  Data.  We  assembled  a 
data  file  of  reported  GPS  bit-hit  anomalies  from 
PC  relational  datab^es  maintained  by  AF  Global 
Weather  Central  (AFGWC)  and  the  National 
Geophysical  Data  Center  (NGDC),  for  the  period 
1  October  1984  to  30  September  1990, 


The  50  OSS/WE  provided  addmonal  GPS 
anomaly  reports  in  hardcopy  for  1  October  1990 
through  31  March  1991;  these  were  manually 
added  to  the  data  file  Table  1  is  an  example  of 
the  GPS  anomaly  infonaation  stored  ui  the  data 
file. 


TABLE  1.  GPS  satellite  anomaly  data.  Latitude  (Lat)  and  longitude  (Lon)  data  is  given  in  decimal 
degrees.  Longitude  is  degrees  east  of  Greenwich.  The  satellite  orbital  inclination  (Inc)  is  rounded  to  the 
nearest  degree. 


Satellite 

ID 

Anomaly 

Date 

Start 

Time 

(Z) 

End 

Time 

(Z) 

Start 

Time 

(M 

Lat 

Lon 

Alt 

(km) 

Inc 

(Deg) 

GPS  5113 

84/10/01 

12:01:00 

12:20 

10:47 

51.0 

341,7 

17600 

63 

GPS  9783 

87/11/04 

10:17:00 

11.00 

19:01 

00.6 

130.2 

20042 

63 

2.2  Space  Environmental  Data.  We  obtained 
particle  fluxes  observed  from  geostationary 
satellites  and  geomagnetic  indices  from  the 
USAFETAC  SESS  Climatic  Database  (SCDB). 
Solar  radio  flux,  x-ray  reports,  and  GOES  particle 
fluxes  were  obtained  from  NGDC.  GOES  particle 


data  was  available  for  the  period  1  January  1986 
through  31  March  1991  only.  All  other  data  types 
were  available  for  1  January  1984  through  31 
March  1991.  Table  2  describes  the  space 
environmental  data  used  in  the  study. 
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TABLE  2.  Spac«  Environmental  Variables.  is  a  dimensionless  index  diat  provides  a  linear 
measure  of  the  level  of  disturbance  of  the  geomagnetic  Geld.  X-ray  flares  were  observed  by  the  GOES 
satellites. 


Environmuntal  Variable 

Data  Type 

Source 

Units 

24-Hour  Aj 

Geomagnetic  Index 

SCDB 

None 

1  3-Hour  a , 

Geomagnetic  Index 

SCDB 

None 

30-300  KeV  Electron  Flux 

Satellite  Particle 

SCDB 

Particles/Steradian/Sec/MeV 

95-300  KeV  Electron  Flux 

Satellite  Particle 

SCDB 

Particles/Steradimi/Sec/MeV 

.2-2  MeV  Electron  Flux 

Satellite  Particle 

SCDB 

Particles/Steradian/Sec/McV 

10.7  Cm  Radio  Flux 

Solar  Flux 

NGDC 

lO'^*'  Watts/M=/Hertz 

X-ray  Flare  Report 

1-S  Angstrom  Flux 

NGDC 

Watts/M' 

>2  MeV  Electron  Flux 

Satellite  Particle 

NGDC 

Particles/S  teradian/Sec/Me  V 

>1  MeV  Proton  Flux 

Satellite  Particle 

NGDC 

Particles/Steradian/Sec/MeV 

>5  MeV  Proton  Flux 

Satellite  Particle 

NGDC 

Paiticles/Steradian/Sec/MeV 

>10  MeV  Proton  Flux 

Satellite  Particle 

NGDC 

Particles/Steradian/Sec/MeV 

>50  MeV  Proton  Flux 

SatelHte  Particle 

NGDC 

Particles/Steradian/Sec/MeV 

2.3  GPS  Orbital  Data.  We  estimated  GPS 
latitudes  and  longitudes  missing  from  the  satellite 
anomaly  reports  using  GPS  orbital  prediction 
models.  Archived  GPS  orbital  variables  were  used 
to  initialize  the  models.  Amospace  Corporation 
provided  the  GPS  orbital  data  and  two  orbital 
prediction  models,  PCSOAP  and  DGEN. 


PCSOAP  is  the  Personal  Computer  Satellite  Orbit 
Analysis  Program;  it  models  the  dynamics  of 
artiGcial  satellites  orbiting  the  earth.  DGEN  is  a 
mainframe  model  that  does  integrations  of  the 
equations  of  motion  to  predict  the  locations  of 
GPS  satellites. 
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3.  METHODOLOGY 


3.1  Datebas«  Creation.  Using  two  databases  (a 
GPS  anomaly  database  and  an  environmental 
database  covering  the  period  1  October  1984 
through  31  March  1991),  the  state  of  the  space 
environment  during  GPS  anomalies  could  be 
compared  to  the  general  state  of  the  space 
environment.  The  GPS  and  environmental  data 
collected  for  this  study  were  put  into  the  two 
datab^es.  The  GPS  anomaly  start  times  were 
matched  to  the  corresponding  environmental  data 


to  form  the  GPS  miomalv  dati^ase,  which  contains 
the  environmental  variables  (mentiored  in  Table 
2)  that  specified  the  space  environment  dunng 
each  GPS  anomaly.  The  time  window  before  the 
start  of  each  anomal}  for  each  environmental 
variable  observation  (as  shown  in  Table  3)  varies 
because  the  amount  of  time  varies  before  different 
parts  of  the  space  environment  (such  as  the 
energetic  particle  or  geomagnetic  environment) 
react  to  solar  activity. 


TABLE  3.  GPS  Environmental  Observations.  The  time  resolutions  described  below  are  the  result 
of  consultations  between  solar  forecasters  at  AFGWC  and  space  analysts  at  USAFETAC;  they  are  the 
variables  used  in  real-time  GPS  anomaly  assessment. 


Environmental  Data  Type 

Observation  Windovr 

X-ray  Flare 

Largest  x-ray  flare  reported  within  24  hours  of  anomaly  start 
time 

24-hour  Ap 

Value  for  calendar  day 

3 -hour  Up 

First  value  reported  after  anomaly  begins 

10.7  cm  Radio  Flux 

Value  for  calendar  day 

Satellite  Particle  Fluxes 

Particle  fluxes  at  the  time  anomaly  begins 

The  environmental  database  is  a  control  database 
separate  from  the  GPS  anomaly  database,  and 
contains  environmental  variables  reported  for  each 
hour  of  each  d^  during  the  period  1  October 
1984  through  31  March  1991.  The  exception  to 


this  is  the  GOES  particle  data,  which  was  not  kept 
in  the  NGDC  archive  prior  to  1  January  1986. 
Table  4  shows  time  resolution  for  the 
environmental  variables  in  the  database. 
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TABLE  4.  EnvIroniDsnUil  Observations.  Tiic  time  resolutions  desenised  bciov.  uert  rci|ucstcd  tn 
the  solar  forecasters  at  AFGWC 


Environmental  Data  Type 

Database  Time  Resolution 

X-rav  Flare 

Laa-gesi  fla/e  m  progress  during  eat’h  hous 

24-hour  A^ 

Value  for  calendar  dsv 

3-hour  a. 

Latest  value  reported  before  start  of  each  hour 

10  7  cm  Radio  Flux 

Value  for  calendar  dav 

Satellite  Particle  Fluxes 

Maximum  5 -minute  value  reported  each  hour 
Minimum  5-minutc  value  reported  each  hour 
.^vcragc  5-mtnoic  value  reported  each  hour 

3.2  Statistical  Comparisons.  Statisucai 
Analysis  System  (SAS)  software  was  used  to  sort 
the  variables  in  the  GPS  anomaly  database  and  the 
environmental  database  into  statistical  distributions 
and  to  plot  each  distnbution  to  allow  differences 
between  the  environment  during  satellite 
anomalies  and  the  normal  space  environment  to  be 
seen. 

3.3  Statistical  Calculations.  SAS  was  also 
used  in  a  multi-step  process  on  the  environmenui 
database  to  quantify  any  correlation  between 
activity  ui  the  space  C’r  ironment  «id  the 
occurrence  of  GPS  anomalies  Since  the  diftenng 
time  resolutions  between  the  GPS  anomaly 
database  and  the  environmental  database  made  it 
impossible  to  calculate  direct  c/'rrelations  between 
the  environmental  variables  and  anomaly 
occurrence,  a  modified  version  of  the 
environmental  database  was  used  for  the  statistical 
calculations.  First,  a  new  element  was  added  to 
the  environmental  database:  whether  or  not  an 
anomaly  occurred  during  each  hour  of  the  penod 
between  1  January  1986  and  31  March  1991 
Even  though  we  have  GPS  anomaly  reports  for 
1984-5,  data  from  these  years  was  not  used  in 
these  calculations  because  GOES  proton  fluxes 
were  not  archived  for  that  period. 


Correlation  coefficicais  were  calculated  between 
ait  the  environmental  vanabtes  and  the 
occurrence  non -occurrence  of  anomalies  during 
each  hour  of  the  studs  penod 

Next,  a  stepw  ise  linear  regression  w  as  done  on  the 
modified  environmental  database  to  detrrmine 
which  environmental  vanablcs  had  the  best 
correlation  to  the  occurrence  of  an  anomaJv 
Since  some  of  the  effects  of  a  solar  disturbance 
take  several  davs  to  reach  the  Earth,  linear 
regressions  were  done  on  the  environmental 
vanablcs  dunng  the  hour  the  anomah  occurred 
and.  24.  48.  and  72  hours  before  the  anomaH 
occurred  In  each  case,  those  environmental 
vanables  selected  b>  SAS  as  the  most  significant 
tn  the  occurrence  of  anomalies  were  noted 

Ftnally.  disenmmant  analysis  was  used  to  estimate 
the  probability  of  anomalies  occumng.  based  on 
the  different  environmental  vfjnables  previousK 
selected  Uiscnroinant  anatysis  produces  a  roatnx 
that  shows  how  well  the  selected  variables  can 
classifv  the  environmental  database  into  two 
categones  either  a  GPS  anomalv  occurred  or  did 
not  occur  Discriminant  analysis  uses  the 
regression  equations  previousK  derived  to  predict 
the  probability  of  an  anomah  occumng  for  each 
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hour  of  the  stmly  period,  and  uses  that  probiU»ilit>  hour  of  the  study  penod,  w>  jmoroalv  occuired 

to  sort  the  study'  data  into  yes/no  categories  for  when  the  probability  was  greater  than  50% 

anomaly  occurrence.  We  assumed  that  for  any 
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4.  LIMITATIONS 


4.1  Percent  of  Anomalies  Report«d.  In  order 
to  do  the  statistical  calculations  described  in 
Section  3.3,  the  assumption  bad  to  be  made  that 
every  GPS  anomaly  that  occurred  between  1 
October  1984  and  31  March  1991  is  in  the  GPS 
anomaly  database.  This  assumption  is 
questionable  for  two  reasons: 

•  There  are  long  periods  of  time,  such  as  the  last 
6  months  of  1988  and  the  first  4  months  of  1989. 
during  which  no  GPS  anomalies  were  reported. 

•  Even  though  the  AFGWC  and  NGDC  satellite 
anomaly  databases  cover  the  same  period  of  time 
and  were  supposed  to  contain  the  same 
information,  they  had  only  five  anomaly  reports  in 
common. 

4.2  Lack  of  Environmental  Data.  The 
magnetosphere  is  a  data-sparse  region;  few 
observations  of  the  space  environment  are  taken 
consistently  there.  For  this  study,  we  used 
whatever  satellite  particle  data  was  available,  even 
though  the  observing  satellites  were  often  in  a 
different  part  of  magnetosphere  than  the  GPS 
satellites  when  an  anomaly  occurred.  The  particle 
data  observed,  therefore,  may  not  have  been 
representative  of  the  space  environment  at  the 
GPS  satellite.  Although  some  limited  datasets  of 
particle  data  observed  by  sensors  on  GPS  satellites 
exist,  USAFETAC  was  not  able  to  obtain  any  GPS 
data  for  this  study. 


4.3  incomplete  Satellite  Reports.  In  using 
orbital  dynamics  models,  missing  ops  locations 
had  to  be  filled  after  the  fact,  this  introducal  the 
possibilit>  of  error  to  the  database  A  more 
significant  problem  in  the  GPS  database  was  the 
lack  of  end  times  for  GPS  anomalies  In  order  to 
complete  this  study,  we  m^e  the  assumption  that 
ail  anomaly  reports  without  end  times  (about  35% 
of  the  database)  were  brief  anomalies  of  less  than 
1  hour  in  duration. 

4.4  Accuracy  of  Anomaly  Timas.  Twenty 
percent  of  the  anomaly  start  times  reported  in  the 
database  were  the  times  that  the  anomalies  were 
reported  to  the  groundstation,  and  not  the  tunes 
that  the  anomalies  actually  occurred;  GPS 
anomalies  with  start  times  exactly  on  the  hour, 
along  with  those  with  start  times  between  2200 
and  2212Z,  fall  into  this  category  All  sudi 
anomalies  were  eliminated  from  the  database  and 
the  statistics  were  recalculated.  Dropping  the 
questionable  anomalies  changed  the  results  of  the 
statistical  analysis  by  less  than  2%  and  had  no 
effect  on  which  environmental  variables  were 
selected  as  significant  predictors  of  anomalies. 
Statistical  comparions  between  the  daily  space 
environment  and  the  occurrence  of  GPS  anomalies 
also  showed  no  difference  in  the  enviroiunental 
varibales  sleeted  as  significant  or  in  the  results  of 
the  analysis. 
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S.  RESULTS 


5.1  Statistical  Comparison  Rosults.  The 
Appendix  contains  charts  of  distributions  of  each 
environmental  variable  in  the  GPS  anomaly 
database,  along  with  the  corresponding  variable 
firom  the  mtvironmental  database.  Comparisons 
between  the  environment  during  anomalies  and  the 
general  environment  show  several  contradictory 
features.  Some  of  the  charts  show  a  higher  than 
expected  frequency  of  occurrence  for  higher 
values  of  the  variable.  The  x-ray  flare  chart,  and 
many  of  the  high-energy  proton  charts,  show  a 
cluster  of  values  in  the  "tail”  of  the  distributions, 
implying  that  for  some  anomalies  the  environment 
is  disturbed.  However,  many  charts  also  show 
higher  than  expected  frequency  of  occurrence  for 
lower  values  of  the  variable.  The  10.7  cm  flux, 
the  dp,  the  24-hr  Ap,  x-ray  flares,  and  some  of  the 
low-energy  electrons  show  this  feature.  These 
charts  show  that  any  relationship  between  the 
space  environment  and  satellite  anomalies  is 
complex  and  not  very  well  understood. 

5.2  Statistical  Correlation  Results.  Linear 
correlations  between  the  occurrence  of  satellite 


anomalies  and  the  environmental  variables  were 
calculated.  The  correlation  coefficients  indicate 
the  amount  of  linear  dependence  between  each 
environmental  variable  and  the  occurrence  of 
satellite  anomalies.  If  the  environment  were 
ideally  correlated  with  the  occurrence  of 
anomalies,  the  correlation  coefficient  would  be  -1 
or  +1,  showing  total  linear  dependence  A 
correlation  coefficient  of  zero  indicates  the 
variable  are  linearly  independent  of  anomaly 
occurrence.  Figure  3  shows  the  correlation 
coefficients  for  the  environmental  database 

All  the  correlation  coefficients  are  extremely  low  , 
showing  hardly  any  linear  dependence  between  the 
occurrence  of  anomalies  and  the  space 
environment.  In  the  case  of  environmoital  data, 
a  correlation  coefiFicient  of  >.6  or  <-.6  would  be  a 
sign  of  a  strong  correlation.  The  best  correlations 
seen  above  are  for  hourly  minimum  >50  MeV 
protons  and  the  10.7  cm  flux.  These 
environmental  vari^les  will  be  frequoitly  selected 
during  the  stepwise  linear  regression  as  being 
statistically  significant. 
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Environmontal  Variabtoa 

Corralation 

Environmantal  Variabiaa 

Corralation 

>50  MeV  Proton  Flux  min 

.088 

10.7  Cm  Radio  Flux 

-087 

>5  MeV  Proton  Flux  avg 

.066 

>10  MeV  Proton  Flux  min 

066 

>5  MeV  Proton  Flux  min 

.066 

>10  MeV  Proton  Flux  avg 

.065 

>5  MeV  Proton  Flux  max 

.064 

>10  MeV  Proton  Flux  max 

.063 

>50  MeV  Proton  Flux  avg 

.063 

>1  MeV  Proton  Flux  avg 

061 

>1  MeV  Proton  Flux  avg 

.059 

>1  MeV  Proton  Flux  max 

.053 

.2-2  MeV  Electron  Flux  min 

.037 

>50  MeV  Proton  Flux  max 

036 

24-Hour  Ap 

-.019 

3-Hour  Sp 

-Oil 

>2  MeV  Electron  Flux  max 

.009 

>2  MeV  Electron  Flux  avg 

009 

>2  MeV  Electron  Flux  min 

.008 

X-ray  Flare  Report 

.007 

.2-2  MeV  Electron  Flux  max 

-.005 

95-300  KeV  Electron  Flux 
min 

004 

30-300  KeV  Election  Flux 

max 

-.003 

30-300  KeV  Electron  Flux 
avg 

-.003 

30-300  KeV  Electron  Flux 
min 

-.003 

95-300  KeV  Electron  Flux 
avg 

-.001 

95-300  KeV  Electron  Flux 

max 

-.001 

.2-2  MeV  Electron  Flux  avg 

-.001 

FIGURE  3.  Pearson  correlation  coefficients  for  environmental  database.  Correlated  variables 
are  listed  from  left  to  right  in  order  of  descending  importance. 


5.3  Discriminant  Analysis  Results.  We  tried 
sevmsl  different  discriminant  analysis  techniques 
to  select  which  environmental  variables  are  most 
signilRcant  in  showing  differences  in  the  space 
environment  between  cases  in  which  GPS 
anomalies  do  or  do  not  occur.  Tables  containing 
classification  matrices  will  show  how  well  certain 
groups  of  environmental  variables  did  in 
classifying  die  environmental  database  into 
categories  of  occurrence/non-occurrence  of 
anomalies.  Skill  scores  were  also  used  to  quantify 
the  accuracy  of  the  environmental  variables  in 
classifying  anomalies. 


5.3.1  Skill  Scores.  A  skill  score  is  a  statistical 
formula  used  to  determine  if  a  categorical  (yes/no) 
weather  forecast  technique  will  work  better  than 
random  chance.  For  diis  study,  die  Hanssen  and 
Kuiper  discriminant  "V"  score  was  used  to 
determine  the  relative  accuracies  of  the 
environmental  variables  used  to  classify'  the 
environmental  database  into  yes/no  cases  of 
anomaly  occurrence.  This  skill  score  is  considered 
the  best  for  rare  events  like  satellite  anomalies. 
The  discriminant  "V"  score  (VDS)  was  calculated 
using  Equation  2. 
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VDS 


AD-BC 


(2) 


where  A,  B,  C,  and  0  as  shown  m  the  mamx; 
Note  that  A,  B,  C,  md  D  are  ohso^ation  couats 
and  not  percentages. 


QPS  ANOMALIES 

Expoetad 

Occumd 

No 

Ya» 

No 

A 

B 

Yes 

C 

VDS  ranges  between  -1  (no  skill;  worse  than 
random  diance)  mid  -^1  (totally  accuntte).  A  VDS 
of  ZNO  meant  the  environmental  variables  are  no 
better  tham  chance  in  detennining  if  an  anomaly 
may  occur.  For  die  purposes  of  this  study,  a 
VDS  higher  dim  .3  indicates  that  the 
envuonmentid  variables  used  to  classify  anomalies 
do  significandy  better  dian  chance. 

•  5.3.2  Ciassifleation  Matricas.  Skill  scores 

alone  are  not  enough  to  judge  die  accuracy  of 
environmental  variables  in  forecasting  or 
classifying  events.  Classification  matrices  ate  also 
used  to  provide  more  detailed  information  on 
important  factors  like  "capability,"  or  the  number 


of  missed  events  (events  not  forecast  thm 
occurred),  and  "reliability"  (the  number  of  "fidse 
almms,"  or  events  forecut  that  did  not  occur). 
Equations  3  and  4  show  the  formulas  used  to 
calculate  reliability  (REL)  and  capidiility  (CAP). 


REL 


D 

(B*D) 


(3) 


CAP 


D 

(C*D) 


where  B,  C,  and  D  are  observation  counts  as 
shown  in  the  matrix  on  the  next  pige. 
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Reliability  shows  the  number  of  conectly 
ctassified  unomaly  events  over  the  total  numbm'  of 
evmts  classified  as  anomalies.  The  higher  the 
reliabili^,  the  fewer  false  alarms.  CtgfobilHy 
shows  die  number  of  correctly  classified  anomalies 
vs  the  total  number  of  miomalies  that  occurred. 


The  higher  the  cap^ility,  the  fewer  misted 
anomalies.  The  regression  equmions  used  to 
classify  all  the  environmoital  observations  into  the 
matrices  in  this  section  ate  listed  in  Figure  1.  The 
mmrices  in  this  report  use  the  format  shown  in 
Table  5. 


TABLE  5.  Claasiflcation  Matrix  Formal  The  vdiues  in  the  matrices  are  numbm 
of  obsen'ations  used  to  calculate  dcill  scores,  reliability,  and  capdiility,  as  well  as  the 
percent^es  of  observations  for  easier  comparison  between  groups  of  variables.  The 
pm:enu^e$  rutming  horiaontaliy  across  the  matrix  (the  expected  yes  and  no)  add  up  to 
100%.  The  skill  scores,  reliability,  and  cigiability  for  alt  groups  of  environmental 
variables  will  be  shown  in  Table  9  at  the  end  of  this  section. 


OPS  ANOMALIES 

Exp«ct»d 

Occumid 

No 

Yos 

No 

#  Correct 

#  False  Alarms 

%  Correct  non-anomaly  evmits 

%  False  Alarms 

Yes 

a  Missed  Anomaly 

#  Correct 

%  Missed  Anomalies 

%  Correct  anomaly  events 

•  5.3.3  Environm«nt  Ourtng  AnomalNis.  The 

results  shown  in  Table  6  are  for  two  groups  of 
eavironmoBtal  variables  for  die  hours  that  GPS 
anomalies  occurred.  Both  groups  had  identical 
"V"  skill  scores  of  .34,  which  border  on  being 


significantly  better  than  chance.  There  is  some 
vtfiation  between  the  groups  in  the  percentages  of 
false  alarms  and  missed  anomd;ie$,  but  not  enough 
to  be  significant. 


TABLE  6.  Classification  of  tha  anvironmant  during  OPS  anomaiias.  The 

variables  listed  below  each  matrix  are  in  their  ordm*  of  statistical  importance  as  selected 
by  the  SAS  stepwise  discriminant  analysis  procedure. 


Group  1 


25420 

15600 

62 

38 

313 

797 

28.3 

71.7 

Qrouo  1  vartablaa:  Hourly  minimum  >S0  MeV  protons;  10.7  cm  flux;  x-ray  flux; 
hottriy  minimum  .29-2  MeV  electrons;  hourly  minimum  95-300  KcV  electrons;  hourly 
minimum  >10  MeV  protons;  hourly  maximum  >50  MeV  protons;  hourly  maximum  >10 
MeV  protons;  hourly  memi  >5  MeV  protons. 


Group  2 


26217 

17718 

59.7 

40.3 

306 

900 

25.4 

74.6 

Group  2  varlabloa:  Hourly  minimum  >50  MeV  protons;  10.7  cm  flux;  x-ray  flux; 
hourly  minimum  >10  MeV  protons;  hourly  maximum  >50  MeV  proU>ns;  houriy  maximum 
>10  MeV  protons;  houriy  mean  >5  MeV  protons. 


In  Group  2  we  tried  to  improve  the  accuracy  of 
the  discriminant  analysis  by  dropping  die  low 
enmgy  electrons,  which  seemed  to  be  less 
significant  dtan  die  high  eaargy  protons,  x-n^ 
flmiM,  and  10.7  cm  radio  flux.  The  Group  2 
results  had  a  slightly  lower  percentage  of  missed 
anomalies  but  a  higher  false  alarm  rate.  These 
results  imply  diat  low  energy  electrons  by 
themselves  do  not  have  a  significant  effect  on  the 
occurrence  of  anomalies.  This  is  in  line  with  the 
results  of  die  stepwise  regression  technique,  which 
iiidicided  Aua  &e  correlation  betwemi  any  one 


variable  and  the  occurrence  or  non-occurrence  of 
anomalies  is  weak. 

5.3.4  The  Environment  24-72  Hours  Before 

GPS  Anomalies.  Since  some  effects  of  solar 
disturbances  can  take  several  days  to  effect  the 
near-Eaith  space  environment,  we  did  statistical 
stedies  of  the  environment  before  anomalies 
occurred  to  determine  if  any  environmental 
variables  had  a  significant  effect  on  anomaly 
occurrence.  The  results  shown  in  Table  7  are  for 
environmental  variables  observed  24,  4S,  and  72 
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hours  t^ore  each  anomaly.  As  the  anomaly  time 
Ijf^s  further  behind  the  mviroiimmital  observation, 
the  variables  selected  chmige  but  the  matrix 
percentiles  and  the  skill  scores  degrade  only 
slightly.  The  "V"  skill  score  was  .3  for  the 


24-hour  case,  .29  for  die  48-hour  ctae,  mid  .28  for 
the  72-hour  case.  In  all  three  cues,  &e  skill 
scores  ud  the  percenti^te  of  correctly  classified 
anomalies  are  lower  than  in  the  cases  at  the  times 
of  die  anomalies  described  in  Section  5.1  3 


TABLE  7.  Classification  of  ths  snvironmsnt  24-72  hours  bsforo  OPS 
anomaliss.  The  three  groups  of  variables  listed  below  are  in  their  order  of  statistical 
importance  as  selected  by  die  SAS  stepwise  discriminant  andysis  procedure.  Since  the 
percentages  in  the  classification  matrices  changed  by  less  than  1  percent  for  the  three 
cases,  one  matrix  with  the  percentages  rounded  off  to  the  nearest  d^ree  is  shown 


Groups  3-5 


22573 

18415 

55 

43 

299 

834 

26 

74 

Group  3  vsrtabiss  f24-hours  prsviousi;  10.7  cm  flux;  hourly  mean  >lMeV  protons; 
244ottr  Ap;  hourly  minimum  .29-2  MeV  electrons;  x-ray  flux;  hourly  maximum  30-300 
KeV  electrons. 

Group  4  vartshlos  f4S4tours  prsviousi;  10.7  cm  flux;  hourly  minimum  >lMeV 
protoiu;  hourly  minimum  .2-2  MeV  electrons;  hourly  mmimum  >5MeV  protons;  hourly 
minimum  >10MeV  protons;  hourly  mean  >lMeV  protons;  3^our  a^. 

Group  5  vsrisbiss  f72-hours  prsviousi:  10.7  cm  flux;  hourly  minimum  .29-2  MeV 
electrons;  hourly  minimum  >S0  MeV  protons;  3^our  a^;  hourly  maximum  >lMeV 
protons;  houriy  minimum  95-300  KeV  electrons. 

Stadsticalty,  the  10.7  cm  radio  flux  is  die  most 
important  environmmital  variable  in  all  three 
cases.  The  others  don't  contribute  as  much,  which 
is  why  die  matrix  percentages  and  skill  scores  ate 
similar  in  all  three  cases.  As  the  time  lag  before 
aiomaly  occurrence  lengthens,  x-ray  flare  reports 
become  less  significant  and  die  geomagnetic  index 
(a^)  becomes  more  significant. 

5.3.5  Environmsnt  During  Solar  Minimum 
and  Maximum.  The  period  of  time  covered  by 
the  statistical  study  (1986  to  1991)  includes  a  solar 
cycle  mbmitm  in  1986  and  a  maximum  in  1989. 

We  divubd  die  dat^ase  into  two  cases:  die  solar 


minimum  years  1986-88  and  the  solar  maxuoum 
years  1989-91  to  see  if  the  general  level  of  solar 
activity  had  any  influence  on  the  correlation 
between  the  environmental  variables  and  the 
occurrence  of  GPS  anomalies.  Table  8  shows  the 
results  for  Group  6  (significrmt  enviroiunental 
variables  during  the  solar  minimum  period)  and 
Group  7  (significant  environmental  variables 
during  the  solar  maximum  period).  The  "V"  skill 
score  for  Group  1  is  slightly  lower  (.30)  dian  for 
the  entire  study  period  (see  Section  5.1.3),  while 
die  score  for  Group  7  is  the  lowest  found  in  this 
study  (.12). 
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TABLE  8.  Clastmcation  of  tho  onvironmont  during  solar  min  and  max.  The 

vari^les  listed  below  each  m^x  are  in  their  order  of  statistical  importance  as  selected 
by  the  SAS  stepwise  discriminant  analysis  proctKlure. 


Group  6 


11843 

11470 

50.8 

49.2 

171 

642 

21.0 

79.0 

Group  8  variablas:  10.7  cm  flux;  hourly  minimum  >50  McV  protons;  hourly 
minimum  .29-2  MeV  electrons;  hourly  minimum  95-300  KeV  electrons,  hourly  maximum 
30-300  KeV  electrons;  hourly  mean  30-300  KeV  electrons;  hourly  minimum  30*300  KeV 
electrons. 


Group  7 


19049 

258 

98.7 

1.3 

290 

43 

87.1 

12.9 

Group  7  vartablea;  Hourly  maximum  >5  MeV  protons;  hourly  maximum  >50  MeV 
protons;  hourly  mean  >50  MeV  protons;  x-ray  flux;  hourly  mean  >1  MeV  protons;  hourly 
mean  >10  MeV  protons;  hourly  maximum  >10  MeV  protons;  hourly  minimum  >10  MeV 
fHotons. 


The  10.7  cm  radio  flux,  50  MeV  protons,  and  low 
energy  electrons  show  the  greatest  influence  as 
anomidy  predictors  in  the  solar  minimum  case.  In 
die  solar  maximum  case,  high  energy  protons  (>1 
MeV  and  up)  and  the  x-n^  flux  are  preferred 
variables. 

Similar  pmcentages  occur  for  the  solar  minimum 
classification  matrix  and  results  for  the  entire 
ptewd  of  die  study  (see  Ttdile  6).  The  percentage 
of  missed  anomalies  dropped  by  7%,  but  at  the 
expemse  of  a  higher  percentage  of  false  alarms. 
However,  the  solar  maximum  matrix  shows  very 
d^(»mit  results.  The  fdse  alarm  rate  during  solar 
max  is  practically  nonexistent,  but  87%  of  the 
anomalies  that  occurred  were  not  correctly 
classified,  resultii^  in  a  much  tower  skill  score. 


This  result  is  in  line  widi  an  empirical  finding  by 
solar  forecasters;  diat  is,  during  periods  of  quiet 
solar  acdvily  it's  easier  to  correlate  the  occurrence 
of  occasional  solar  disturbances  with  particular 
satellite  anomalies.  During  periods  of  high  solar 
activity,  there  are  so  many  disturbances  occurring 
almost  simultaneously  that  it's  difiTicult  to 
determine  which  (if  any)  might  be  the  cause  of  a 
particular  satellite  anomaly.  The  higher  levels  of 
energetic  particles  saturate  the  satellite  sensors  and 
affect  how  satellites  react  to  changes  in  the 
environment. 

Table  9  shows  die  "V"  skill  scores,  capability,  and 
reliability  for  all  seven  groups  of  environmental 
variables  studied. 
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TABLE  9.  Summary  of  results  for  all  groups  of  environmental  variables. 


1  Group 

"V  Skill 
Score 

Capability 

Score 

Reliability  I 
Score  1 

1 ;  9  Variables 

.34 

.72 

.05 

2;  7  Variid>les 

34 

.75 

.05 

3;  24-hrs  Previous 

.30 

.74 

.04 

4:  48-hrs  Previous 

.29 

.74 

.04 

1  5:  72>hrs  Previous 

.28 

.73 

,04 

6;  Solar  Minimum 

.30 

.79 

.05 

7;  Solar  Maximum 

.12 

.13 

.14 

Groups  1  and  2,  die  general  enviroomental  groups, 
had  the  best  skill  scores.  The  others,  with  the 
exception  of  Group  7,  had  scores  only  slightly 
worse.  Group  7,  the  solar  maximum  group,  had  a 
significantly  lower  skill  score  than  the  other 
groups,  due  to  the  large  number  of  missed 
anomalies.  The  c^ability  scores  show  that  when 
anomalies  occur,  all  groups  except  for  Group  7  did 
well  in  correctly  classifying  those  events.  For 
Group  7,  large  numbers  of  missed  anomalies  made 
die  cqiability  score  low.  Reliability  scores  are 
very  low  for  all  groups,  because  the  number  of 
false  alarms  are  alwrtys  very  high  compared  to  the 
number  of  correctly  classified  anomalies.  The 
large  difference  in  size  (two  orders  of  magnitude) 
between  the  number  of  observations  in  which  no 
anomaly  was  reported  and  the  number  in  which 
an  anomaly  did  occur,  biases  the  results. 


Group  7  did  somewhat  better  in  reliability  only 
because  the  number  of  false  alarms  was  much 
lower  (only  one  order  of  magnitude  larger)  than 
for  the  other  groups. 

5.3.6  Geomagnutle  Latitude  and 
Occurrence  of  Anomalies.  In  addition  to  the 
normal  environmental  variables,  we  also  examined 
whether  or  not  a  satellite's  position  had  any 
bearing  on  the  occurrence  of  an  anomaly.  We 
calculated  hourly  positions  for  GPS  satellites  5113, 
5114,  and  5118  for  the  years  1987-91.  Figure  6  is 
a  comparison  between  distributions  of  corrected 
geomagnetic  latitude  during  anomalies  versus  the 
distribution  of  corrected  geomagnetic  latitudes 
along  the  satellites'  orbit. 
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Figure  4.  Distribution  of  Positions  for  GPS  5113,  5114,  and  5118  (corrected  geomagnetic  latitude). 


Certain  geomagnetic  latitudes  have  a  higher 
percentage  of  anomalies  than  would  be  expected 
from  the  normal  orbital  track.  One  area  of  high 
anomaly  frequency  is  in  the  low  latitudes  of  the 
Northern  Hemisphere,  possibly  due  to  the 
equatorial  anomaly.  Another  area  of  relatively 
frequent  anomalies  is  in  the  Northern  Hemisphere 
polar  region;  this  may  be  caused  by  the  orientation 
of  the  magnetic  field  lines  into  the  Earth. 
However,  when  statistical  analysis  was  done  on 
the  effect  of  the  environment  on  one  satellite  (GPS 
51 13),  geomagnetic  latitude  was  not  selected  as  a 
statistically  significant  predictor  of  satellite 


anomalies.  The  skill  score  and  classification  matrix 
were  virtually  identical  to  that  for  the  general 
environment  (Groups  I  and  2  in  Section  5.3.3). 
Examination  of  Figure  6  raises  other  questions 
The  high  frequency  of  anomalies  seen  near  the 
equator  and  pole  in  the  Northern  Hemisphere  is 
not  seen  in  the  Southern  Hemisphere.  Are  the 
satellites  commanded  more  in  the  Northern 
Hemisphere,  triggering  the  occurrence  of  an 
anomaly  (such  as  a  sudden  discharge)?  The 
extremely  small  size  of  the  GPS  5113  anomaly 
database  used  in  analysis  (241  anomalies)  may 
have  also  handicapped  our  analysis. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  Conclusions.  The  objective  of  this  study 
was  to  examine  a  group  of  GPS  anomalies 
collected  over  a  period  of  time  and  determine  if 
there  were  relationships  between  commonly 
observed  space  environmental  variables  and  the 
occurrence  of  anomalies. 

We  found  the  relationship  between  environment 
and  anomaly  occurrence  to  be  somewhat  better 
than  chance,  but  could  find  no  definite  link 
between  any  one  environmental  variable  and  the 
occurrence  of  anomalies.  The  Group  2  variables 
did  best  in  categorizing  anomalies  in  general,  but 
the  large  variability  seen  in  environmental  data 
over  the  solar  cycle  and  the  very  low  correlations 
between  environmental  variables  and  satellite 
anomalies  make  it  risky  to  use  these  results  for 
more  than  general  guidance.  Unless  these  results 
are  tested  on  independent  data  and  found  to  be 
valid,  they  should  not  be  used  to  evaluate 
individual  anomalies. 


6.2  Recommendations.  Better  documentation 
of  GPS  anomalies  is  needed.  More  consistent 
reporting  of  positions,  start  and  end  times,  and 
whether  or  not  the  satellite  was  commanded  when 
an  anomaly  occurred,  might  help  separate 
environmental  factors  from  engineering  factors  in 
classifying  anomaly  causes.  The  organizations 
involved  in  satellite  operations  need  to  make  sure 
that  all  anomalies  are  reported  and  stored  in  the 
anomaly  database  now  maintained  by  the  AF 
Space  Forecast  Center.  The  amount  and  quaiitv  of 
space  environmental  data  will  improve  over  the 
next  10  years  as  more  sensors  and  better  models 
come  into  operational  use.  Once  the  space 
environment  at  the  location  of  each  satellite  can  be 
specified,  a  follow-on  to  this  study  several  years 
from  now  might  yield  more  definitive  results.  The 
importance  of  space  systems  necessitates 
continuing  study  of  space  climate  and  its  influence 
on  system  effects. 
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APPENDIX 


Solar  and  Geomagnetic  Comparison  Charts 

A-i  Distribution  of  Solar  Flares,  X-rj^  Levels. 

A-2  Distribution  of  Daily  10.7  cm  Radio  Flux  10-22  Watts/Sq  M/Hz 

A-3  Hourly  Distribution  of  GOES  Protons,  >  50  MeV 

A-4  Hourly  Distribution  of  GOES  Protons,  >  5  MeV 

A-5  Maximum  Hourly  Distribution  of  GOES  Protons,  >10  MeV 

A-6  Minimum  Hourly  Distribution  of  GOES  Protons,  >10  MeV 

A-7  Hourly  Distribution  of  GOES  Protons,  >50  MeV 

A-8  Distribution  of  3-Hour  Ap  Geomagnetic  Index 

A-9  Distribution  of  24-Hour  Ap  Geomagnetic  Index 

A- 10  Hourly  Distribution  of  Energetic  Electrons  .29  -  2  MeV 

A-1 1  Hourly  Distribution  of  Energetic  Electrons  95  -  300  MeV 

A-12  Daily  Satellite  Anomalies  1986 

A-13  Daily  Environmental  Data  1986 
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A-l.  Distribation  of  Solar  Flares.  X-ray  Lesds. 
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A-2.  Distribution  of  Daily  10.7cm  Radio  Fins  10-22  Watts/Sq  MfHz. 
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A-3.  Hourly  Distributioa  of  GOES  Protons,  >  50  MeV, 
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A-4.  Hourly  Distribution  of  GOES  Protons,  >  5  MeV. 
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A-«.  Minimum  Hourly  Distribution  of  GOES  Protons,  >  10  McV 
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A-7.  HoHriy  Distribution  of  GOES  Protons,  >  50  MeV 


0  10  20  30  40  50  60  70  80  90  100110120130>140 


■  Anomaly  S  3-Hour  Ap 


A-8.  Distribution  of  3-Hour  Ap,  Geomagnetic  Index 
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A-9,  DistribHtioa  of  24-HoMr  Ap,  Geomagnetic  Index. 
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A-IO.  Hourly  Diitribution  of  Energetic  Electronr  .29-2  MeV. 
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A-ll.  Honriy  Distribatioa  of  Enersetic  Electrons,  95-300  KeV. 
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A-12.  Daily  Satellite  Anomalies  1986. 
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COMNAVOCEANCOM.  Code  N312,  SMnnie  Space  Or.  MS  39529-5000  .  2 

NAVOCEANO  (Rue^  Rueeum),  Bidg  8100,  Rm  203D,  Stennia  Space  Or,  MS  39522-5001  .  2 

NAVOCEANO,  Code  9220  (Tony  Ortolano},  Stennie  Space  Ctr,  MS  39529-5001  .  1 

Mauiy  Oceanographic  Ubrary,  Naval  Oceanography  OfBee,  Stennia  Space  Ctr,  MS  39522-5001  .  1 

Naval  Raeearch  Laboratory,  Montaray,  CA  93943-5006  . . .  1 

Naval  Reaeareh  Laboratory,  Code  4323,  Washington,  DC  20375  .  1 

Naval  Reaeareh  Laboratory  (Dr  Rley),  Code  4180,  Washington,  DC  20375  .  1 

Naval  Postgraduate  School,  Chmn,  Dept  of  Meteorology.  Code  63,  Monterey,  CA  93943-5000  .  1 

Naval  Air  Warfare  Centor-Weapone  Division,  Geophysical  Scisnces  Branch,  Code  3254,  Point  Mugu,  CA  93042-5001  .  1 

Chief,  APG  Met  Team,  Bldg  1134,  Attn:  AMSTE-TC-AM  CAB,  Aberdeen  Proving  Ground.  MO  21005-5001  .  1 

Atnwsphertc  Sciences  Laboratory  (SLCAS-AS-1 3  10-2c),  White  Sands  Missile  Range,  NM  88002-5501  .  1 

TECOM  Atnxts  Sci  Div,  AMSTE-TC-AA  (MacBiain),  While  Sands  Missae  Range,  NM  88002-5504  .  1 

VWHte  Sands  Met  Team,  AMSTE-TC-AM  (WS),  White  Sands  Miss8e  Range,  NM  88002-5501  .  1 

USA  TECOM,  ATTN;  AMSTE-TC-AM  (RE)  TECOM  Met  Team,  Redstone  Arsenal.  AL  35898-8052  .  1 

Director,  U.S.A.-CETEC,  Attn;  GL-AE  (VMittmarsh),  Fort  Behroir,  VA  22060-5546  .  1 

•  Technical  Ubrary,  Ougway  Proving  Ground,  Dugway,  UT  84022-5000  .  1 

NOAAMASC  Ubrary  MC5. 325  Broadway,  Bouider,  CO  803036328  .  2 

NOAA  Ltorary-EOC4W5C4,  Attn  ACQ,  6009  Executive  Btvd,  RockviBe,  MO  20852  .  1 

NOAA/NESDIS  (Attn;  Nancy  Everson,  EmA22),  World  Wsather  Bldg,  Rm  703,  Washington,  DC  20233  .  1 

NGOC,  NOAA,  Mall  Code  EA3C4, 325  Broadway,  Boulder,  CO  80333-3328  .  1 

NIST  Pubs  Production,  Rm  A635,  Admin  Bldg,  GaRhersburg,  MO  20899  .  1 

NASA-MSFC-ES44,  Attn;  Dale  Johnson,  Huntsvile.  AL  35812-5000  .  1 

NASA-MSFC-ES44,  Attn:  Gwenevere  Jasper,  HuntsvINe,  AL  358126000  .  1 

Cmc-FDAC,  Cameron  Station,  Alexandria,  VA  22304-6145  .  2 

AUULSE .  Maxwell  AFB.  AL  361 1 2-5564  .  1 

AWSTL,  Scott  AFB,  IL  822256438  .  35 
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